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f*s • We consider the problem whether a nonparametric zero-curvature 

representation can be embedded into a one-parameter family within 
the same Lie algebra. After introducing a computable cohomological 
obstruction, a method using the recursion operator to incorporate the 
parameter is discussed. 
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Ph ■ 1 Introduction 



Despite of four decades' development of soliton theory, a general algorithm 
to recognize integrability still remains elusive. By an integrable system in 
two independent variables x, y we shall mean one that can be reproduced as 
the compatibility condition 

m 

O; Ay-B x + [A,B}=0 (1) 

■^J- ■ of the auxiliary linear system 

o 

W x = AV, m y = BV. 



([47, 38]). Matrices A, B are supposed to belong to a matrix Lie algebra g. 
The q- valued 1-form Adx+B dy is called a zero curvature representation. To 



5_i ■ be indicative of integrability, matrices A, B must depend on what is called 



the spectral parameter. 

Sometimes nonlinear systems come equipped with a nonparametric zero- 
curvature representation. In geometry of immersed surfaces, the Gauss- 
MainardL-Codazzi equations are always of the form (1) while the linear 
problem (2) is supplied by the Gauss- Weingarten equations [46]. Integrable 
systems of this kind relevant to physics were first noticed by Lund and 
Regge [26]. For numerous classical and new examples see [4, 38]. For a dif- 
ferent source (the nonparametric Lax pairs [22]) see [24]. 

Thus we face the "spectral parameter problem" or whether a given non- 
parametric zero curvature representation can be embedded in a non-trivial 



one-parameter family. An arbitrary number of removable parameters can be 
always introduced through the gauge transformation 

A' = SxS' 1 + SAS~\ B' = SyS- 1 + SBS' 1 , 

where S is an invertible matrix belonging to the Lie group G associated with 
the Lie algebra g. The gauge transformation is a left action of the group 
of G- valued functions on the set of all solutions of system (1). However, 
parameters introduced that way play no role in soliton theory. 

In the geometric context, the majority of effective approaches link the 
spectral parameter with a deformation parameter (Bobenko [4]). A nonre- 
movable parameter can be also often identified with the group parameter 
of a suitable point symmetry, see, e.g., [25, 41]. The early nineties [23, 5] 
witnessed a development of algorithmic procedures based on existing al- 
gorithms to compute symmetries. To detect integrable cases, Cieslihski [5] 
compares dimensions of the point symmetry algebras of the input system 
and of its covering [20] induced by the zero curvature representation. If 
these are inequal, a candidate for integrability is discovered. Obviously, un- 
detected remain the cases when the spectral parameter cannot be identified 
with a group parameter. The first counterexample to be investigated thor- 
oughly was that of inhomogeneous nonlinear Schrodinger system. This led 
Cieslihski to consideration of extended symmetries [6, 7] that operate within 
a class of equations. This method is suitable for classification problems, but 
results can depend on the initial class of equations considered (the extended 
symmetry must operate within this class) . In this paper we present a general 
non-existence criterion for the spectral parameter problem within a given Lie 
algebra, a given jet order, and a given order of power expansion. Our main 
result can be viewed as a "restricting" procedure which says when a nonre- 
movable parameter cannot be inserted. We also consider a new "extending" 
procedure that seeks the way to insert the parameter, thus trying to fill the 
space left by the restricting procedure. The spectral parameter problem is 
solved (in certain jet order) when the results of restricting and extending 
procedures meet. 

2 Characteristic elements 

In this section we recall few necessary facts from the geometric theory of 
partial differential equations [21] and zero curvature representations [27]. 
Consider an arbitrary system of partial differential equations 

E l (x,y,u i ,ul,ul...)=0 (3) 



of an arbitrary order. By Uj we shall denote the derivative of u specified by 
a symmetric multiindex I in x, y. In the jet space with coordinates t, x, vlj 
we consider total derivatives 

Equations (3) and their differential consequences DjE = determine a 
submanifold (called diffiety), which we shall denote by £. Here Dj denotes 
the composition of the operators of total derivative D x , D y (in any order, 
since they commute). Finally, operators D x , D y admit restriction to the 
diffiety. 

A zero curvature representation is defined to be a 1-form a = Adx+B dy 
such that 

D y A - D X B + [A, B] = (4) 

on £ . Two zero curvature representations Adx + B dy and A' dx + B' dy are 
said to be gauge equivalent if one can be obtained from the other by the 
gauge transformation 

A' = D X S S- 1 + SAS-\ B' = D y S S' 1 + SBS' 1 . 

The gauge equivalence explains why system (4) is so much underdetermined 
(has twice as many unknowns as equations), which is the main obstacle to 
successful search for its solutions. To compute zero curvature representa- 
tions, we can choose between the classical Wahlquist-Estabrook method of 
prolongation algebras [45, 9] (see [14] for implementation), its extensions 
such as [10], the method of coverings by Krasil'shchik and Vinogradov [20], 
the remarkable development by Igonin [17], the Sakovich method of cyclic 
bases [39, 40], and also the method of characteristic elements [27, 28]. Being 
a modification of the well-known characteristic function method to compute 
conservation laws, the latter approach tries to fix the gauge and copes with 
classification problems reasonably well. 

Given a zero curvature representation Adx + B dy of system (3), the 
characteristic elements can be computed using the Sakovich formula [39]. 
Namely, we necessarily have D y A — D X B + [A, B] = ^jxCfDjE 1 with 
suitable matrices C[ , and then the Ith characteristic element Cm is given 
by 

C {l) = ^(-D)jCf . (5) 



Here D x = D x — [A, ■ ], D y = D y — [B, ■ ], and Dj is their obvious composi- 
tion (in any order since D x ,D y necessarily commute once A, B form a zero 
curvature representation). Finally, \z denotes the restriction to £. 

Example 1. The mKdV equation ut + u xxx — 6u 2 u x = has an st(2)-valued 
zero curvature representation Adx + B dt with 

_ (u A\ (~u xx + 2v? - AXu 2Xu x + 2Xu 2 - 4A 2 



1 -uj ' ~ V ~ 2u x + 2n 2 - 4A u xx - 2u 3 + 4An 
Actually, we have Dt(A) — D X (B) + [A, B] = (ut + u xxx — 6u 2 u x )-Cn\, where 



C ^ \0 -1 

By formula (5), Cn\ is the characteristic element. The only nonzero sum- 
mand in (5) corresponds to the empty multiindex J. This is typical for zero 
curvature representations computable by the Wahlquist-Estabrook proce- 
dure. 

Under gauge transformation, characteristic elements transform by conju- 
gation, which allows us to bring one of them into the Jordan normal form J. 
This usually leaves some gauge freedom, which can be exploited to reduce 
one of the matrices Ai by gauge transformation with respect to the relatively 
small stabilizer subgroup Sj of J . This very much resembles the well-known 
"wild" problem of normal forms of pairs of matrices with respect to simul- 
taneous conjugation. In reality our problem is simpler, since zero curvature 
representations of use in soliton theory never take values in a resolvable sub- 
algebra, which allows us to cut branches that fall into such a subalgebra. 
In the case of Lie algebra s[(2) and s((3) all the normal forms were found 
in [28] and [43], respectively. For convenience we tabulate them in the two 
columns headed k = of Table 1 and 2. The remaining columns, headed 
k > 0, will be explained in Section 3. 

The general case of Lie algebra s((n) still remains largely unexplored 
except when the Jordan normal form J consists of a single Jordan block, 
which has been completely investigated by Sebestyen [44]. 

All possible zero curvature representations (4) along with their charac- 
teristic elements (5) now can be found from the determining system 



D y A-D x B + [A,B]\ e = 0, 



£(-£) 



i,i 



dE 



du 



a 



fe°(0 



(6) 
0, v ; 



£ 



Table 1: Normal forms in the case of algebra sl(2) 
C A 



k = 


k>0 


k = 


k>0 


(o-") 


(o 4) 


fan ai 2 \ 
V 1 -an) 


fan ai 2 \ 
V -an) 



1\ /0 0\ / a i2 \ / 012 
J l c J 1 021 J I a 2 i 



derived in [27]. With A,B,Cm restricted to their respective normal forms, 
the determining system (6) enjoys the following properties: 

- is a system of differential equations in total derivatives; 

- has the same number of unknowns as equations, akin to systems de- 
termining symmetries and conservation laws; 

- is quasilinear in A, B, and linear in Cm. 

Solution of system (6) is essentially algorithmic given an upper bound on 
the jet order of the unknowns. Because of nonlinearity in A and B, solution 
can be troublesome even with the help of computer algebra. 

Remark 1. Wahlquist-Estabrook type zero curvature representations have 
the order equal to the order of the equation minus one. The order of their 
normal form can be twice the order of the equation minus one at worst. 

Let us remark that the procedure does not need to fix the gauge com- 
pletely. The general solution of system (6) may still depend on a removable 
"false" parameter. The parameter originates from linearity of system (6) in 
the characteristic elements as demonstarted in the following example. 

Example 2. Consider the case of the Lie algebra sl(2) and the Jordan 
normal form 

of the characteristic element. The 1-parameter gauge group 




s" 



S=L „-i , s 7^0, (7) 



preserves the normal form of A and multiplies J by s 2 . By linearity, if 
A, B, Cm is a solution of system (6), then A, B, s 2 Cm is a solution as well. 



Table 2: Normal forms in the case of algebra sl(3) 

C A 

k = k > k = k > 

'an ai2 Oi 3 \ /on ai2 ai 3 N 

1 022 023 J I a 2 2 023 

y Q31 1 Q33/ \«31 «33 y 

1 \ / 

021 022 1 J 021 a 2 2 

V«31 «32 033/ \031 032 33 , 









'on 012 \ /an 012 

021 a 2 2 1 a 2 i 022 



y0 3 i ^32 O33/ \03i 032 CJ33, 

ai2 \ / ai2 ai3 N 

021 a 2 2 1 J a 2 i 022 

\a 3 i a 33 / \a 3 i a 32 a 33/ 






012 

021 ^22 

\a3i a 33/ 
012 on 
021 

\031 a 32 0; 





'an a 12 1 \ / ai 3 

021 a 22 J a 2 i 022 a 2 3 

V«31 a 32 O33/ \03l 032 033 / 

ai2 \ / 012 ai3 N 

021 022 023 J I 021 ^22 ^23 

^3i a 33 / \a 3 i a 33 , 



All matrices are supposed to be traceless (e.g., c 3 = — ci — c 2 ). 



If Cm = J and A is normalized according to Table 1, then A, B,s 2 Cn\ 
is an unnormalized solution, but SAS" 1 , SBS^ 1 , Ss 2 Cn\S~ 1 is a normal- 
ized solution. Thus (7) introduces one false parameter into the solution of 
system (6) under normalization. 

So far, the only complete classification (with unbounded jet order) ob- 
tained is [30] (second order evolution equations possessing a s((2)-valued zero 
curvature representations; none admitting an essential parameter, hence all 
nonintegrable as naturally expected). 

3 Restriction procedure. Formal removability 

In the previous work [29], a nontrivial horizontal gauge cohomology group 
H 1 was associated with zero curvature representations depending on a non- 
removable parameter. Here we elaborate the idea further in terms of power 
series in the spectral parameter. 

Given a nonparametric zero curvature representation Aq,Bq, consider all 
possible expansions 

oo oo 

A(X) = J2 AiX, B(\) = Y, BiW (8) 

i=0 j=0 

around zero. Obviously, ^o = A(0), Bq = B(0). 

Upon inserting expansions (8) into the zero curvature condition (4) we 
obtain 

D y A k - D x B k + Y,[A U Bj] = (9) 

i+j=k 

for all k > 0. For k = this equation says just that Aq, Bq is a zero curvature 
representation. For k = 1, equation (9) reads 

D y A 1 - D x B l + [A u B ] + [Aq, B ± ] = (10) 

(according to [11], this equation also plays an important role in the geometry 
of surfaces immersed in Lie algebras). Denoting D x = D x — [Aq,-], D y = 
D y — [Bq, ■ ], equation (10) can be rewritten as 

D y A 1 = D x B 1 (11) 

and interpreted cohomologically [29]. 



Recall that the horizontal gauge complex [29] consists of modules g® A fe 
of horizontal forms with coefficients in the Lie algebra g. Here horizontal 
essentially means that the forms contain differentials dx, dy only and A is 
simply the additive group of functions. The operator d acts on an arbitrary 
g- valued function F € g <8> A as 

dF = D x Fdx + D y Fdy G g® A 1 , 

while for an arbitrary g- valued horizontal 1-form U dx + V dy £ g <8> A 1 we 
have 

d(U dx + V dy) = (D X V - D y U) dxAdy £g® A 2 . 

It is easy to see that dod = once Aq, Bq is a zero curvature representation. 
The group 

1 Keid : g® A 1 — >g<g)A 2 



Imd : g® A — > g <£> A 1 

is called the first horizontal gauge cohomology group of the zero curvature 
representation Aq,Bq. Elements of Kerc? are called cocycles while those of 
Imd are coboundaries. 

We show that H 1 = implies removability of the parameter in a strict 
sense. To circumvent the issue of convergence, we adopt the following defi- 
nition. 

Definition 1. The parameter A is said to be formally removable if every 
summand in the expansions (8) except Aq,Bq can be locally annihilated by 
a gauge transformation with respect to some matrix function in A. 

The following proposition is easy to prove. 

Proposition 1. Let A(\),B(\) be a zero curvature representation depend- 
ing analytically on a parameter A. Let the first horizontal gauge cohomology 
group H 1 with respect to Aq, Bq be zero. Then A is formally removable. 

Proof. Equation (11) says that A\ dx + B\ dy is a cocycle, hence also a 
coboundary if H 1 = 0. This means that A\ dx + B\ dy = dS\ for a suitable 
g- valued matrix function Si, i.e., 

A 1 = D X S 1 -[A ,S 1 ], 

B\ = DySi — [Bq, Si}. 



Let S(X) = E + S\\. The gauge action with respect to S{\) will keep the 
terms Aq,Bq, but remove the terms Ai,B±. This is readily checked using 
the expansion (E + Si A)" 1 = E - S\X + SfX 2 - SfX 3 +- ■ ■ , which has a 
positive radius of convergence in every open set where E + S\ has bounded 
eigenvalues, see Gantmacher [13, Ch. V, § 4, Thm. 2]. This means that 
A±,Bi can be locally annihilated. Otherwise said, if H 1 = 0, then without 
loss of generality we can assume that A\ = B\ = 0. 

Now we can proceed by induction. If Ai = Bi = for all i < k, then 
equation (9) with k > 2 reads 

D y A k - D x B k + [A k , B ] + [A , B k ] = 0, 

meaning that A k dx + B k dy is a cocycle, hence a coboundary, i.e., equal to 
dS\ for suitable S k . Reasoning as above, one easily sees that A k ,B k can be 
locally gauged out through the gauge matrix E + S k X k . This concludes the 
proof. □ 

Corollary 1. A zero curvature representation Aq,Bq such that the first 
horizontal gauge cohomology group H 1 is zero cannot be a member of an an- 
alytic one-parameter family that would depend on a formally nonremovable 
parameter. 

A method to compute the cohomology group H 1 of a given zero curvature 
representation Aq,Bq follows from [29, Prop. 7], which we reproduce below 
in a slightly simplified form: 

Proposition 2 ([29]). Given a zero curvature representation Aq dx + B^ dy, 
the 1 -form A\ dx + B\ dy is a cocycle if and only if matrices 

constitute a zero curvature representation: 

D y AW - D X BW + [AW , £W] = o. (13) 

Moreover, two cocycles are cohomological (differ by a coboundary) if and only 
if the corresponding zero curvature representations (12) are gauge equivalent 
with respect to a gauge matrix of the block triangular form 

sln - { E s I) < M > 

with the unit matrix E at the diagonal positions. 



It follows that cohomology classes of cocycles Ax dx + Bx dy can be de- 
termined by solving equation (13) modulo gauge equivalence with respect 
to (14). If A\ = 0, B\ = for all normalized solutions to (13), then H 1 = 0. 
However, there can also be a kind of "false" solutions to equation (13). 

Example 3. Continuing Example 2, consider the gauge matrix 



S 



1 + ciA + c 2 A 2 +- 


2 







(l + ciA + c 2 A 2 +---)~\ 

'I +f]A + c 2 A z +--- 

1 - ciA + (c? - c 2 )A 2 +- 

with constant coefficients Cj. In Example 2 we already demonstrated that 
S introduces a false parameter not detected by the characteristic element 



method. What is then the action of S on the matrix Aq 
By expanding SAqS" 1 in the powers of A we get 

SAaS' 1 = A + 



an 

021 



ai2 

-an 







-2ciai 2 V 
2cia 2i / 



+ 







{c\ + 2c 2 )a 2 i 


2cia 2 i 



(3c? -2c 2 )a i2 \ ., 

o |A+ ' 



-2ciOl 2 





! . If this Ai exhausts solutions to (13), 

v v ^ciu, 2 i u / 

again. 

= 0, then Proposition 1 applies. Otherwise, A\,B\ are candidates 



In particular, Ai 
then H 1 -- 

UH 1 

for the first coefficients of the Taylor expansions (8). How to check whether 
they are good? 

It is easily verified that the system of equations (9) indexed by k = 
0, . . . , m is satisfied if and only if the block triangular matrices 



(15) 





(A ••• 0\ 




(Bo 







°) 


^H _ 


A 1 A '•• ': 

: '•• ■■. 


, B^ = 


Bx 


Bo 









Um ••• A x A J 




\"m 




Bx 


bJ 


constitute a zero curvature representation, i.e., 










D V AW - 


-D x BH + [iH )B [ 


m ]] = 0. 











(16) 

Hence, if Ak,Bk are already known for all k < m, then (16) determines 
whether expansions (8) can be extended one step further. 
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Proposition 3. If, for some m, no solution A^ m ',B^- m ' of equation (16) 
exists, then there is no possibility to extend the expansions (8) beyond the 
first m terms. 

Propositions 1 and 3 reduce the spectral parameter problem to that of 
computation of zero curvature representations modulo a gauge group. The 
substantial benefit is that the problem becomes linear, being reducible to 
solution of a system (6) written on matrices A*- 1 ', B*- 1 ' given by (12), resp. 
on matrices A'-" 1 ', B' m ' given by (15). The system (6) is then indeed linear in 
all unknowns including Ak, Bk, k > 1 (being still nonlinear in Aq, Bq, which 
are supposed to be known, however). 

To solve system (13) or (16) by the characteristic element method, we 
need to construct the respective normal forms. The characteristic elements 
of the zero curvature representation A<- k >, B^ k ' are block triangular matrices 
of the same form 



/C(j)c 



C 



[k] 







c, 



(Oi u (0o 



C7, 



\C(i)k 



c, 



(1)1 



\ 



C(i)oJ 



(17) 



Every diagonal block Cn\ coincides with the corresponding known charac- 
teristic element of the input zero curvature representation Aq, Bq. The gauge 
transformations considered should be with respect to matrices of the form 



s [k] 



(E 
Si E 



\S k 



o\ 







Si E) 



In the case of the Lie algebras s((2) and sl(3) the normal forms for C^ and 
Ak, k > 0, are tabelated in Table 1 and 2 above. 

Being a determined system of equations in total derivatives, system (6) 
is routinely solvable if one sets an upper bound on the jet order of the un- 
knowns (see the end of Section 2). It seems natural to assume that insertion 
of the paramater does not require increasing the jet order of the zero cur- 
vature representation. However, this is not necessarily the same thing as 
setting the order of Aj~, B^ equal to that of Aq, Bq, since normalization can 
increase the jet order (see Remark 1 and the example in the next section). 



11 



4 An example. The inhomogeneous nonlinear Schro- 
dinger equation 

In this section we consider Cieslihski's counterexample [1, 6, 7, 8], the inho- 
mogeneous nonlinear Schrodinger equation 

Qt = (fq)xx + 2gr, 

Pt = -(fp)xx-2pr, (18) 

r x = f(pq)x + zfxpq, 

where f(t,x) is an arbitrary function. Compared to Cieslihski's formula- 
tion [6, 7, 8], we have put p = q and made the formal change of coordinates 
t «-» it, preserving integrability. The system is known to be integrable if 
/ = fi(t)x + fo(t) is linear in x; otherwise it is believed to be nonintegrable 
(see op. cit. and references therein). 

The initial nonparametric first order sl(2)-valued zero curvature repre- 
sentation is D t Ao — D x Bq + [Aq, Bq] = with 

valid for all functions f(t,x). 

Let us demonstrate how the procedure outlined in the previous section 
detects the expected integrable cases in two steps. To start with, we find the 
characteristic element of Aq dx + Bq dy to be the triple of s((2)-matrices 

1\ „ / 0\ „ (-\ 



c «°- Vo oj' c w°-{-i oj' C ^° v 1 

Indeed, the expression D t Ag — D x Bq + [Aq, Bq] is equal to the sum of C(j\q 
multiplied by jth equation (18) each, similarly as in Example 1 above. 

To proceed further, we construct the 4x4 matrices (12) and (17), 
k = 1, from the already known blocks Aq,Bq, C^q, CV 2 )oj C(3)o an d the 
blocks A\,Bi, C(i)i, C(2)i> C(3)i Y e * to be found. These matrices are subject 
to system (6). At this point there are essentially two natural ways to pro- 
ceed further depending on which pair of the unknown matrices we choose to 
normalize. 

Choice 1. According to the second row of Table 1 above, matrices CmQ 
and Aq are already in their respective normal forms. However, Cmo is a 
single Jordan block, hence a "false" parameter should be expected (see Ex- 
ample 2). 

12 



According to Table 1 above, normal forms for Cr\\\ and A\ are 



a 



(i)i 




c 



Al 



a 12 
a.21 



Then there is no room for further normalization and matrices B\ , C { 
have to be left arbitrary. Therefore, the 4x4 matrices 



(2)1,C(3)1 





^° 


9 C 


^ 








A® = 


< 


c 

112 C 




q 


; 






V«21 


-i 


0/ 






/ r 


(/<?): 


e o 





BW = 


(fp)x 

Oil 


— r 

012 



r 



(fq)x 




\ 021 


"fell 


(fp) 


x -r 




/O 1 


0\ 






G (i) " 






\c 



1 
0/ 


1 






I ° 








0\ 


G (2)- 


-1 

C(2)ll 




C(2)] 




.2 






; 




\ C (2)21 


-C(2)] 


LI "I 


°/ 




/-l 








0\ 


G (3)- 




C(3)ll 


1 

C(3)] 




.2 "I 










\ C (3)21 


-C( 3) ] 


Ll 





V 





involve 12 unknown functions, namely on,ai2, on, 012, 621, c, C( 2 )n, C( 2 )i2, C( 2 )2i, 
c (3)iii c (3)i2i c (3)2i- The determining system (6) consists of the same number 
12 of independent linear equations in total derivatives. Assuming A\,B\ of 

we obtain a nontrivial solution 



the second order at least (cf. Remark 1 


A 1 =\ p 2 ), B 1 = 

Vo 0/ 


P 



(fp), 



p 



+ 2/0 



ifph 

p 



(20) 
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(and its constant multiples) of system (6), valid without any constraint on 
f(t,x). We also get a solution of the form shown in Example 3, which cor- 
responds to a false parameter. 

This implies that H l is always nonzero. This means that all inhomoge- 
neous nonlinear Schrodinger equations admit a nonremovable parameter up 
to the first order. 

To see whether the power expansion can be prolonged further, as a next 
step we consider the 6x6 matrices A™, B^ 2 \ C,\, CLl, C,l from the already 
known blocks Ai,Bi, Cnv, C(2)j, C^, i < 1 and the yet unknown blocks 
A 2 , B 2 , C(!) 2 , C( 2 )2, C( 3 ) 2 . Normalizing A 2 , C( 3 ) 2 in the same way as A 1 ,C^- )1 
above, we consider the corresponding system (6). However, this new system 
turns out to be incompatible unless 

d 2 f 

I?- - <21) 

By Proposition 3, condition (21) is necessary for the initial zero curvature 
representation (19) to admit a nonremovable parameter, assuming the jet 
order at most two. 

Choice 2. According to the first row of Table 1, Ct^Q is in a normal 
form, whilst Aq is not. Normalizing ^o increases the jet order of Aq,Bq to 
two. This leads us immediately to computations on the second order jet 
level, while avoiding the "false" parameter solution. 

5 Extension procedures. Recursion operators 

The procedure allows us to compute every term of the Taylor expansion of 
A(X),B(X). In the case of polynomial dependence on the spectral parame- 
ter we get a closed-form result after a finite number of steps. Otherwise we 
are left with a truncated Taylor expansion, which is, however, not entirely 
useless. In principle, one can apply the characteristic element method to 
compute the full zero curvature representation ^4(A),.B(A). Solving the non- 
linear system (6) by triangularization [16] usually involves heavy branching, 
but knowing the truncated Taylor expansion (8) permits cutting off many 
branches immediately. 

The symmetry method [5] and its extension [6, 7, 8] by Ciesliiiski take 
advantage of linearity of the problem. Remarkably enough, certain aspects 
of this method admit a cohomological interpretation, which will be discussed 
elsewhere. 
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Apart from zero curvature representations, two-dimensional integrable 
systems usually possess infinite hierarchies of symmetries, generated by re- 
cursion operators. Given a zero curvature representation or a Lax pair, the 
recursion operator 1Z can be derived by various methods, see [12, 40] and 
references therein. Another relation was observed to exist between the zero 
curvature representation and the inverse (1Z + Aid) -1 [3, 31, 33], which even 
extends to certain multidimensional systems [34]. This, however, requires 
interpretation of the recursion operator as a Backlund autotransformation 
of the linearized equation (Papachristou [37], Guthrie [15]). 

Although this interpretation first appeared in multidimensional integra- 
bility theory [37], it turned up to be extremely fruitful in two-dimensional 
case as well. Guthrie's motivating question was how to apply a recursion 
operator to a symmetry in a safe way. The dangerous point is that the com- 
mon pscudo differential form of a recursion operator [36] involves a purely 
formal inverse D~ 1 = D' 1 . The usual answer is to interpret H = D~ 1 F 
as F = D X H, but this can easily lead to errors ('bogus' symmetries [15]). 
Guthrie's idea was to provide also the value G = DfH, where t is the other 
coordinate. This suffices to determine the pseudopotentials up to an integra- 
tion constant. The integration constant then can be safely omitted since it 
only adds 9t(0) to the result. Sergyeyev [42] discussed another interpretation 
problem absent under Guthrie's approach. 

As is well known, infinitesimal symmetries of the system (3) can be put 
in the vertical form 

i,i 1 

where / runs over all multiindices over t, x, coefficients U l being functions 
defined on 6 . The vector field (22) is a symmetry if and only if the functions 
U l satisfy the system 



1,1 1 



--:£ E i(U)\ £ . (23) 



Formally, (3) and (23) can be written together as a system of partial differ- 
ential equations 

E l = 0, 

Eg*-* 

i,i 1 
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where the components IP of a symmetry are additional unknowns. Sys- 
tem (24) will be called the linearized system. This allows us to interpret 
recursion operators as Backlund autotransformations for solutions U l of the 
linearized system (24) (Papachristou [37]). 

Example 4. Considering the linearized Burgers equation 

Ut = u xx + uu x , 

U t = U xx + uU x + u x ll, 

it can be easily checked that the well-known ^-dependent recursion operator 

U' = W x + \ {tu + x)U + \ (1 + tu x )W, 

W X = U, 

W t = U x + uU, 

is actually a Backlund autotransformation. 

General recursion operators of system (3) are to be sought among Back- 
lund autotransformations of the linearized system (24), with pseudopoten- 
tials subject to first-order linear systems of the form (2). The coefficient ma- 
trices A, B are supposed to depend on u] and Uj now. Hence the pair A, B 
constitutes a zero curvature representation over the linearized system (24). 

Numerous examples [2, 3, 31, 33, 34] (as well as unpublished ones) sup- 
port the following conjecture: 

Conjecture 1. For every recursion operator 9\ of a system (3) there exists 
a finite- dimensional Lie algebra g and a zero curvature representation A, B 
of the system (3) such that the pseudpotentials of 9\ can be assembled in a 
single Q-valued nonlocal variable $ subject to conditions 

<f> x = [A,$]+e A (U), <f> t = [B,<f>]+e B (U). (25) 

Here U denotes a seed symmetry, [-, -] the commutator in g, and £ the 
componentwise Frechet derivative 



k,I 



du) 



For 'conventional' recursion operators the algebra g is typically solvable 
and mostly abelian, in which case the zero curvature representation A, B 
reduces to a collection of conservation laws (possibly nonlocal) [32]. E.g., 
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in Example 4 the algebra is 1-dimensional, ut = (u x + 2 u2 ) x being the 
corresponding conservation law. 

On the other hand, the inverse (SH + Aid) -1 of a conventional recursion 
operator D\ + A Id is usually associated with the A-dependent zero curvature 
representation A, B of system (3) (however, it may happen that the algebra 
is solvable). This is another expression of the conventional wisdom that 
a recursion operator 1Z yields a zero curvature representation through the 
eigenvalue problem 9\U = XU. 

What we have obtained is a way from a A-independent zero curvature 
representation to a recursion operator TZ, the inverse (JZ + Aid) -1 , and its 
associated A-dependent zero curvature representation. To make it into a 
working algorithm we must further restrict the form of the recursion oper- 
ator. 

Conjecture 2. For every integrable system (3) there exists a zero order 
recursion operator with diagonal matrix, i.e., of the form 

m(U) = (W) = (cjtT + a)* j ) (26) 

with *$> as in (25), where c^,o*- are functions on £. 

The nonlocalities ^ can take values in an extension of the algebra Q 
of the initial zero-curvature representation. Leaving this aspect aside, we 
propose the following procedure. 

Procedure 1. Given an initial zero curvature representation A = Aq, B = 
B of (3), 

1. write down the linearized system (24) and the associated covering (25); 

2. find all operators (26) such that U' = 9\(U) satisfies the linearized 
system (24); 

3. solve equation U' = 9\(U) + XU for U in terms of U' to obtain the 
inverse recursion operator (SH + Aid) -1 ; 

4. transform (9^ + Aid) -1 to its Guthrie form and identify the associated 
zero curvature representation (which occurs in adjoint representation here). 

This procedure was powerful enough to lead to all recursion operators 
published in the works [2, 3, 31, 33, 34]. However, not always a recursion 
operator could be associated with every value of the spectral parameter (for 
an instance see [31]). Consequently, for certain A-independent zero curvature 
representations the above procedure fails even though the corresponding A- 
dependent family exists. This phenomenon must be further investigated. 
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Nevertheless, in the following section we show that our procedure copes 
well with the zero curvature representation of the inhomogeneous nonlinear 
Schrodinger equation. 

6 The inhomogeneous nonlinear Schrodinger equa- 
tion continued 

Continuing the example of the inhomogeneous nonlinear Schrodinger equa- 
tion, we consider the nonparametric sl(2)-valued zero curvature representa- 
tion (19) valid for all functions f(t,x). 
Step 1. The linearized system (24) is 

Qt = (fq)xx + 2qr, 
Pt = ~(fp)xx -2pr, 
Tx = f(pq)x + 2f x pq, 
Qt = (fQ)xx + 2(rQ + qR), 
Pt = -(fP)xx-2(rP+pR), 
Rx = f{qP + pQ)x + 2f x {qP + pQ). 

The covering (25) is determined by equations 

o v q\ fipu 

-2q ^12 
-2p 0/ W 

(fp) x -(/?) 

(fq) x -It 

-(/p)x 2r 

where ip %3 are components of an s((2)-matrix 

Vll ^\2 





(27) 



Step 2. The equations to be solved are 

D t Q' - fD xx Q' - 2f x D x Q> - (f xx + 2r)Q' - 2qR' = 0, 

D t P' + fD xx P' + 2f x D x P' + (f xx + 2r)P' + 2pR' = 0, (28) 

D X R' - fqD x P> - (fq x + 2f x q)P> - fpD x Q' - (fp x + 2f x p)Q' = 
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where, according to (26), we assume 

p> = C pP + opVn + dp ipi2 + a 2 pip2i, 
Q' = cqQ + aqtpn + aqipn + Oq^i, 



R' = c R R + a^Vn + 4Vi2 + 4^2 



Functions cp,cq,cr, ap,OQ,o^ are the unknowns and can depend on t, x, 

p,q,r and their derivatives p x ,q x ,r t , p xx ,q xx ,r tt , System (28) can be 

solved routinely. As a result we obtain that %K(U) is a constant multiple of 
U unless f xx = again. On the other hand, if / = fi(t)x + fo(t), then we 
obtain a nontrivial solution 



P' = gi P + ^21, 

Q' = giQ + i>i2, 

R' = giR + fiipu + fpipi2 + /#2i, 
where 



(29) 



5i = / fi(r)dr. 



(30) 



The recursion operator given by formulas (29) and (27) can be applied to 
an arbitrary seed symmetry. However, generated symmetries will be nonlocal 
as a rule, which is usually the case with inverse recursion operators. 

Step 3. To invert this operator, we reinterpret formulas (29) and (27) 
so that P', Q' , R' are given and P, Q, R are to be found. From equation (29) 
we get 



P 



P'-iP 



21 



.91 



Q 



Q'-V 



12 



9i 



,R 



R' ~ (/lV'll + fptpl2 + /#2l) 



.91 



(31) 



while ipij will satisfy a system of the same form (25). The term Aid is 
absorbed in the integration constant of (30). Herefrom we can reconstruct 
the corresponding zero curvature representation as 



A{\) 



251 



V 



-p 



1 



B(X) 



25? 

(/P)* + - 

51 



(/«)* 



_/9\ 

51 

25?/ 



(32) 



19 



where g\ is the integral (30) and the integration constant for g\ serves as 
the spectral parameter A. This result is gauge equivalent to that obtained 
by Ciesliiiski [6, 7]. 

That the parameter is nonremovable can be seen from the expansion 
of A in powers of A. Indeed, A\ dx + B\ dy is then gauge equivalent to the 
previously computed generator (20) of H 1 . In particular, the upper and lower 
bounds established in Sections 4 and 6 coincide, meaning that the answer 
is complete (within the Lie group sl(2) and the second jet order of the zero 
curvature representation) . 

Remarkably enough, the inverse of the recursion operator (31) still fails 
to produce a local hierarchy. 
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